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^ — ' Abstract 

f->i ■ The mechanism of the chiral symmetry breaking is investigated in the strong-coupling Abelian gauge 

I ' theories at finite temperature. The Schwinger-Dyson equation in Landau gauge is employed in the 

C^. real time formalism and is solved numerically within the framework of the instantaneous exchange 

approximation including the effect of the hard thermal loop for the photon propagator. It is found 
that the chiral symmetry is broken below the critical temperature T for sufficiently large coupling a. 
The chiral phase transition is found to be of the 2nd order and the phase diagram on the T — a plane 
is obtained. It is investigated how the structure of the chiral phase transition is affected by the hard 
thermal loops in the photon propagator. 
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1 Introduction 

The chiral phase transition in gauge field theories at finite temperature is an interesting phenomenon in 
many respects. In particular it plays an important role when we deal with the early stage of the universe. 
It has been known for a long time that Abelian gauge theories at vanishing temperature subject to the 
chiral symmetry breaking in the strong coupling region and hence the massless fermions acquire a mass in 
a dynamical way (ij ||, y, ^. This phenomenon has been clarified mostly by using the Schwinger-Dyson 
equation ||]. 

It is quite natural from the point of view of the early universe to introduce the temperature effect in the 
analysis of the above phenomenon and see whether the broken chiral symmetry at vanishing temperature 
is restored at sufficiently high temperature and examine whether the phase transition is of the first order 
or second order. There have been several works dealing with this problem in quantum chromodynamics. 
There have, however, been not many works 0, m, || which have studied this question by taking into account 
of the photon mass coming from the hard thermal loop. Thus we try to present the thorough analysis of 
the chiral phase transition in Abelian gauge theories at finite temperature with due consideration on the 
hard thermal loop by using the Schwinger-Dyson equation. 

In our analysis we deal with the Schwinger-Dyson equation in the real time formalism to introduce the 
temperature P, M, O, n2] . To solve the Schwinger-Dyson equation we confine ourselves to the instantaneous 
exchange approximation. The approximation is found to be valid in the high temperature region. In the case 
of the vanishing temperature the vacuum polarization effect is negligible. For high temperature the vacuum 
polarization becomes enhanced and so within our approximation we are forced to take into account of the 
vacuum polarization function in the photon propagator while dealing with the Schwinger-Dyson equation. 

We make mostly numerical calculations in our analysis of the finite-temperature Schwinger-Dyson equa- 
tion. We find that our numerical solutions are quite stable and the results are proven to be reliable. We 
find the clear signal of the second order chiral phase transition as temperature varies. We investigate the 
thermal photon mass effect in the chiral phase transition. We obtain critical curves on the T ~ a plane with 
T the temperature and a the fine structure constant. 

2 Schwinger-Dyson equation at finite Temperature 

Throughout the paper we work in Abelian gauge theories, in particular quantum electrodynamics, with 
massless fermions. The main purpose of our work is to see whether the broken chiral symmetry for large a 
at vanishing temperature is restored at high temperature and, if so, what is the nature of the chiral phase 
transition, i. e. whether the transition is of the 1st order or the 2nd order. Here a is the fine structure 
constant a = e^ /Air with e the electric charge of fermions. We rely on the Schwinger-Dyson equation in 
the real time formalism of the finite-temperature field theory. In the imaginary time approach the integral 
equation is cast into the infinite-component simultaneous equation and we do not prefer to cut off the 
summation in our calculation. 

We start with the brief summary on the Schwinger-Dyson equation for vanishing temperature T = 0. 
The Schwinger-Dyson equation for the fermion self-energy part S(p) reads within the ladder approximation 
in Landau gauge 

nv) = -^e" J ^r^Siq)r^DlT{p - q), (1) 

where D*J'^^{p — q) is the photon propagator at tree level and the self-energy part S(p) is defined through 

% 1 

'^^^^ " 7/ - S(p) + te " A(p2y_B(p2)+je' (^) 

with S{p) the full propagator for massless fermions, and A{p'^) and B{p^) the invariant functions of p^ 
respectively. By the use of the invariant functions the self-energy part is represented such that S(p) = 
(1 — A^p^))-^ + B{p'^). It has been known for a long time ||l|, g, ^ that Eq. (^ allows a non-vanishing 



solution for B{p^) with A{p'^) = 1 if a > ac = 7r/3. Thus the chiral symmetry is broken for the unusuaUy 
large electromagnetic coupling and the massless fermion acquires the dynamical mass for a > 7r/3. 

At finite temperature in the real time formalism 0, 11^, O, Ol the full propagator for fermions may be 
written in the following form: 



iS{p) = 



(3) 



For simplicity we assume that 

^o(po,b1)=A(po,|p1)-l, ImB(po,b1) = 0. (4) 

It will be seen that the above simplifying assumption does not lead to any inconsistencies. We then have 



iS{p) = 



pf-Mipo,\p\)+ie 
where the mass function M(poj \p\) is defined by 

trReS(p) 



Mipo,\p\) 



trl 



In the closed time path method |l3| the spinor self-energy part is given by the following matrix form: 



^l:^\p) = V-\p,p) (^ '^^P^ _^^^.^^^ ) V-\(3,p), 



(5) 



(6) 



(7) 



where V((3,p) is the unitary matrix which connects the thermal vacuum to the zero-temperature vacuum 
and is given by 



vw,p) 



cos(^ — e(po)sin(y9 
e{po)smip cos(^ 



with 



Accordingly wc find 



and hence 



C0S(/3 ■ 



simp 



exp(-/3|po|/2) 
v/l + exp(-/3|po|)' ~"'^ v/l + exp(-/3|po|)' 

Rel]"(p) = ReI](p), 

trReI]"(p) 



M{paAp\) 



trl 



(8) 



(9) 



(10) 



(11) 



The Schwinger-Dyson equation in the real time formalism is written down in the matrix form within 
the framework of the closed time path method. If the vertex part is approximated by the tree form, the 
1-1 matrix element of the equation is represented only by the 1-1 component of the fermion and photon 
propagator respectively and takes the following form, 



E^^(p) 



(2^ 



r^S'\q)Y^Dllip-q). 



(12) 



Here the spinor two-point function S^^{p) is given by evaluating the 1-1 matrix element of the expression 



tS''''{p)^V{(3,p) 



S{p) 
S*ip) 



ViP,p), 



(13) 



and reads 



iS'\p) = {j^ + M(p,,,\p\)) 



„2 ,,r2/ l^^^.„„^,/3bo| 



+ tt5{p^ ~ M\pQ,\p\))ta.nh'- 



p^-M^po,\p\) '" '"""''' 2 

The photon two-point function D}^],(p) is also given by evaluating the 1-1 matrix element of the matrix 



(14) 



iDfiq) = U{(3,q) (^ ^^'^'^ ^,^^^^ ) UW,q), 



(15) 



where 



cosh^^ , ' ; sinh^^ oxp(-/3|go|/2) _ ^^^^ 

,/l-exp(-/3|go|) v^l - exp(-/3|go|) 

The general form of the photon propagator at finite temperature is well-known and is given by [|ll| 

^" ^^'~ q^-nT{q)+i£ ^^"^ ^ g2 _ n^ (g) + ,^ ^M- * ^2 + ^^ ^2 ' l-^^J 

where Pj^^ and P^^^ are the transverse and longitudinal projection operators respectively and acF is the 
gauge fixing parameter. Below we take the Landau gauge acF — which is consistent with the Ward- 
Takahashi identity within the ladder approximation at T = 0. We then find that 

iDll (p) = ReD^, (g) coth ^ + zlmD^, (?) , (19) 



where 



ReD^,{q) = Rei?^(q)Pj;+Rei5^(q)P^^,, (20) 

ImD^^iq) = ImZ?^(g)Pj,+Imi?^(<7)Pi;. (21) 



with 



ReD^iq) = V -Inin^('?) 



ReP>^(g) = r 



ImD^iq) = V 



lmD^{q) = V 



iq^-UT{q)){q2-nT*{q)) 

neiq' - RelFiqW - Ren^(g))<5((g2 - U^{q))iq^ - n^*(g))), (22) 

-lmn^(g) 
(g2-ni(q))(q2_ni*(q)) 
ns{q' - Ren^(g))(g2 - Ren\q))S{{q^ - n\q)){q^ - n^*(g))), (23) 

q2 -Ren^(q) 

(g2-n^(g))(g2_n^*(^)) 

ns{q' - RelFiq))lmn^{q)d{iq' - ^^ {q)){q^ - n^*(g))), (24) 

q'^ -Rcn^(q) 
(g2-ni(q))(q2_ni*(q)) 
7re(g2 _ Ren^(g))Imn^(g),5((g2 _ n\q)){q^ - n^*{q))). (25) 



Finally the Schwinger-Dyson equation at finite temperature in the real time formalism takes the following 
form 



M{poM) 



1 



-trRe 



-^^ I -i^,r^S^\q)^^lDll{p~q). 



{2nY 



-e" I ^M{qoAg]) 



(27r)4 



V- 



M^qo, \q]) 



(2ReL»^(p -q)+ ReD^{p - q)) coth 



/5bo -90 



+ {2ImD'^{p -q)+ ImD^{p - q))S{q'^ - M^{qo, \q])) tanh 
where it should be noted that the following formulae have been employed, 



P\qo\ 



(26) 



(27) 



3 Instantaneous exchange approximation 

In order to derive informations on the phase structure of quantum electrodynamics at finite temperature 
as much as possible we try to solve the Schwinger-Dyson equation regarding it as an integral equation for 
the mass function M(pQ,\p\). To solve the equation it is inevitable to make some approximation which 
may not give any serious influence on the resulting physical predictions. Throughout the paper we apply 
the instantaneous exchange approximation in which the po dependence of the relevant Green functions is 
assumed to be weak and is neglected. For photons the approximation implies that the vacuum polarization 
function is po independent. As a functional form of the vacuum polarization function we adopt the one 
suggested by the hard thermal 1-loop calculation: 



n^ (90,9)190^0 ^ 0, 



n^ (90, 9)1,0-0 - 2Nml^ = je^T^, 



(28) 



where N is the number of fermion flavors. With this approximation the Schwinger-Dyson equation is 
rewritten as 



Mipo,\p\) 



M{qo,\q\)[2r 



(2^ 
x6iq^-Miqo,\q\f)teinh 



1 



{p-cf}^ 
P\qo\ 



V 



1 



(P - 9)^ + 27Vmp?,2 



(29) 



Since the right-hand side of Eq.(B9) has nopo dependence in this approximation, the mass function M{po, \p\) 
is independent of po and so is written as Af(|p|). Hence we obtain 



Mip) 



a 

2^ 



q'^dq 



M{q) 



2pq 



2 In 



vV + M(9F 

(p + q? 



{p-qf 



In 



{p + qf + iNfUph" 
{p - qY + 2Nmph' 



tanh 



P^ 



M{qY 



(30) 



where we set p = \p\ and q = \q\. In the following sections we analyze Eq. (^) numerically by using 
computers and study the behavior of the mass functions to derive informations on the phase transitions. 



4 Numerical solutions 



We would like to solve Eq. (^0|) by the use of the numerical method. There are several different methods 
available for solving the integral equation (pOf) numerically. Among those methods there exist two standard 
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Figure 1: Typical behavior of the solution of Scliwinger-Dyson equation (|30|) for A^ = 0, a = 1, and p = 0.1 A. 



methods. The one is to discretize the integral in Eq. (pfl) and regard the resulting equation as simultaneous 
equations for the mass function. The other of the two is to start with the suitable trial function for the 
solution and to iterate it until we reach a stable solution. It seems that the latter method is much easier to 
handle and is useful as long as the convergence of the iteration is guaranteed. We would like to adopt this 
latter method in the following arguments. 

We begin with the simplest possible choice for the trial mass function in the iteration. Thus the trial 
function is chosen to be a constant independent of p, 



M{p) = constant. 



(31) 



At each step of the iteration the integration is performed by the use of the Monte Carlo method and 
the integral is cut off at mass scale A. After the first iteration the resulting mass function acquires a p 
dependence and is substituted for the mass function in the integral on the right hand side of Eq. (|30|). 
Repeating this procedure we may have a possible stable result. Whether we have a stable result or not 
should be always checked after sufhciently many iterations. In each calculation in the following we confirm 
the stability of the solution, that is, we obtain the same result starting from the different trial functions. 

We first consider the case without the thermal photon mass A^ = 0. Note here that parameter N plays 
a role of switching on (iV = 1) and off {N — 0) the photon mass and also represents the number of fermion 
flavors. In Fig. 1 the mass function normalized by the cut-off A of the p integration is presented as a 
function of the number of the iterations for the case with p — A/100 and T = O.OIA, 0.20A, 0.25A. Here we 
adopt the unit system with k = 1 where k is the Boltzmann constant. The fluctuations observed in Fig. 1 
are given rise to by the errors in the Monte Carlo integration. 

As is seen in Fig. 1, the resulting mass function becomes stable after about 100 iterations. We then 
push forward our analysis by taking more points for the values of the momentum and obtain the momentum 
dependence of the solution for the mass function. In Fig. 2 the mass function normalized by the cut-off is 
presented as a function of the momentum p nornralized by the cut-off. The mass functions were obtained 
after 1200 iterations. It should be noted here that the convergence of the iterations becomes slower if we sit 
near the critical value of the coupling constant a and temperature T. In this case we need more iterations 
to obtain a stable result for the mass function. 

The case with the thermal photon mass A^ = 1 is studied essentially in the same way as above. We 



M(p) 




Figure 2: Typical shape of the mass function M{p) for A^ = and a = I. 



may also study the case with three fermion flavors N 
present the numerical results on the mass function. 



3 essentially in the same manner. Here we skip to 



5 Chiral phase transitions 

We first wish to observe the behavior of the mass function at some fixed value of p as a function of parameters 
a and T. The a dependence of the mass function M{p) with p = O.IA is shown in Fig. 3 for various fixed 
values of T. Note here that the errors coining from the fiuctuations observed in Fig. 1 are smaller than the 
size of the mark shown for each sample point in Fig. 3. 

As is clearly seen in Fig. 3, the chiral phase transition is of the second order since the fermion mass is 
generated at a critical value of the coupling constant a without any discontinuity. By picking up the critical 
value of a where the mass generation occurs one can obtain the critical parameters ac and Tc respectively. 
If Tc is plotted versus ac, we obtain the phase diagram. We shall come back to this subject later in Section 
6. It should be noted here that the value of the critical coupling constant ac for T = is found to be 
0.44477 ±0.00053 (See Appendix A). This value differs from the well-known value tt/3. The reason for this 
is simple, i. e., we applied the instantaneous exchange approximation in our approach which is not the 
good approximation in the low temperature region. Thus we have to be careful in applying our method and 
consider that our formulation is suitable rather in the high temperature region. 

We can also observe the behavior of the mass function M{p) at some fixed value of p as a function of 
temperature T. The T dependence of M{p) with p = O.IA is given in Fig. 4 for fixed a. Here again we 
clearly observe the second order chiral phase transition in Fig. 4. 

If the thermal photon mass is included {N = 1), the essential feature of the behavior of the mass function 
is more or less the same as in the case without the thermal photon mass. We present the a and T dependence 
of the mass function M{p) are presented in Figs. 5 and 6. 

We recognize that the behavior of the mass function in the low temperature region as seen in Figs. 4 and 

6 is significantly different from the one in the case without the thermal photon mass. The similar behavior 
is observed in the case of iV = 3 although the result is not shown in figures. 
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Figure 3: Dynamical fcrmion mass as a function of the coupling constant a for A'^ = at p = O.lA 



M{p) 



3.5 

3 
2.5 

2 
1.5 

1 
0.5 





1 

t\ A A A A 


1 

A 


A 


1 

A 


A 


1 


:sc X X X X 


X 


X 


X 


"^^^^^^^^ 


_ 


[ ni n n n n n n 


° [Hun 


'^^ 


%^ 




"^x 

X 

x_ 
a = 0.8 O 




1 ' 




1 


n 
n 
n 


a = 1.0 + - 
a = 2.0 n 
a = 3.0 X - 

a = 4.0 A 

1 


.^^^^ 



0.2 



0.4 0.6 

T/A 



0.8 



Figure 4: Dynamical fcrmion mass as a function of the temperature T for A^ = at p = O.lA 
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Figure 5: Dynamical fcrmion mass as a function of the coupling constant a for N — 1 at p — O.lA 
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Figure 6: Dynamical fcrmion mass as a function of the temperature T for A^ = 1 at p = O.lA 




Figure 7: Critical curves for N = 0,N=1,N — 3 and N = oo at p = O.IA. 



6 Critical curve and critical exponents 



By observing the figures obtained in the last section we can directly draw the critical curve on the T — a 
plane. In fact, for the case with A^ = we pick out the values of a where the mass function vanishes and 
plot those values as a function of temperature. The resulting curve is the critical curve for the case with 
N = 0. In order to perform this procedure in a more systematic way and to obtain critical exponents 
simultaneously we apply the following method. 

We fit the function M(p) with n data points Mi, ai and Ti with i = 1, • ■ • n near the critical point by 
assuming the functional form of the function M{p) such that 



M{p)^e^^{a^a,y, 



(32) 



for fixed temperature T and 



M(p)-e^°(^,-^)^ 



(33) 



for fixed coupling constant a respectively. Here CT,Ca,acTTc,i^ and rj are adjustable parameters with ac 
and Tc corresponding to the critical coupling constant and critical temperature respectively and v and rj 
designate the critical exponent. To estimate the values of CT,Ca,ac, Tc, v and 77 we adopt the familiar 
least-squares fit. We choose n data points near the critical point and minimize the quantities 



^[InM, - {v \n{Tc - T,) + Ct)] 



(34) 



with temperature T fixed and 



^[InM, - (r7ln(a, - a,) + C„)]^ 



(35) 



with coupling constant a fixed. 
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Figure 8: Critical exponents v and rj estimated at p — 0.1 A. 



The critical curves for N — 0, N— 1, N~3 and N ~ oo are shown in Fig. 7. Taking the N ^ oo limit 
in Eq. ( |30| ) we find that the dynamical fermion mass for A^ = oo is obtained by replacing a by (2/3)a in the 
case N = 0. Thus the critical coupling constant Oc for A^ = oo is given by ac for TV = multiplied by 2/3. 
The longitudinal mode of the photon propagator acquires the non- vanishing thermal mass for N ^ 0. Since 
the thermal photon mass rriph tends to suppress the gauge interaction, it has an effect of protecting the 
chiral symmetry. As is seen in Fig. 7 the effect of the thermal photon mass enhances the critical coupling 
ac and becomes stronger as T increases. 

In Fig. 8 the critical exponents v and ry are shown as a function of the temperature and coupling constant 
respectively. As is seen in Fig. 8 the critical exponents do not depend on the number of fermion flavors 
A^ guaranteeing the accuracy of our numerical analysis. Taking the average of the critical exponents v for 
a = 1, 2, 3, 4 and r] for T/A = 0.2, 0.4, 0.6, 0.8, 1.0 we obtain 



and 



0.46 for A^ = 0, 
0.50 for A^ = 1, 
0.45 for A^ = 3, 



0.49 for A^ = 0, 
0.51 for N = 1, 
0.49 for AT = 3, 



(36) 



(37) 



respectively. Therefore our results are consistent with v — 1/2 and ry = 1/2 which are in accord with the 
critical exponents in the four- fermion theory ||Tq|. 

In Fig. 9 the A^ dependence of the mass function is presented. Since the longitudinal mode for photon 
propagator acquires the thermal mass for N ^ 0, the photon propagator has only one massless mode. On 
the other hand the photon propagator has three massless modes for A^ = 0. As is seen in Fig. 9 the 
dynamical fermion mass strongly depends on the number of massless modes in the photon propagator near 
the critical temperature where the instantaneous exchange approximation is valid. 

Note that the behaviors of the critical curve and the dynamical fermion mass for finite A^ seem different 
from the results in Ref. 0|. In Ref. ^ the thermal photon mass is introduced both for the longitudinal 
and transversal mode. In this case there is no massless mode for the photon propagator. 
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7 Conclusion 

We have investigated the mechanism of the chiral symmetry breaking in the strong-couphng Abehan gauge 
theories at finite temperature. The Schwinger-Dyson equation in Landau gauge is solved numerically within 
the framework of the instantaneous exchange approximation which is valid at high temperature. The effect 
of the hard thermal loop for the photon propagator on the phase transition has been also clarified. Here 
are our results in order: 

1. The chiral phase transition is found to be of the 2nd order in the high temperature region. Thus the 
physical mass obeys the scaling law of the mean-field type. The effect of temperature on the chiral 
phase transition seems to be equivalent to the fermion loop effect in strong coupling QED at T = 0. 

2. The critical temperature grows linearly as the critical coupling constant grows in the high temperature 
region. This result is opposed to that of Ref. where they claimed that the chiral symmetry is always 
restored at finite temperature no matter how large the coupling constant is taken. The origin of this 
discrepancy lies in the vacuum polarization function. We have taken its transverse and longitudinal 
parts separately (see Eq. (Pq)) while in Ref. [Q they were assumed to be the same. Note that our 
choice of the function gives no constraints on the number of fermion flavors, whereas in strong coupling 
QED at r = there might be the critical value Nc above which there is the chiral symmetric phase 
only ^. 

3. The thermal photon mass reduces the effect of electromagnetic interaction. In the limit of the infinite 
number of fermion flavors the behavior of the mass function in the case with the thermal photon mass 
(A'' = oo) is the same as in the case without the thermal photon mass (A^ — 0) with the coupling 
constant a replaced by 2a/3. 

The next step that we should put forward is to extend our analysis to the one in the low temperature 
region. Along this line there are some difficulties to be overcome: 

1. The vacuum polarization function at low temperature has to be estimated with some suitable approx- 
imation or in numerical calculations. 
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2. Proper approximations have to be found out in order to solve the Schwinger-Dyson equation at low 
temperature otherwise the equation is much complicated. 

We expect that the structure of the chiral phase transition is the same as in the case we have investigated. 
This will be a subject in the forthcoming paper. 
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Appendix 

Here we argue the critical value of coupling constant ac at vanishing temperature within the instantaneous 
exchange approximation, though this approximation is valid only at high temperature. 
We begin with the linearized Schwinger-Dyson equation 



Mip) 



3a 

Anp 



dqM{q) In 



{p + q? 
[p - qf 



(A.l) 



which is taken from Eq.(pO) by setting T = and neglecting the mass function in the denominator. If we 
take only the first term in the expansion of the logarithm. 



In 
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2n+l 



ip -qf ^ 2n -t- 1 
it is easily found that the Schwinger-Dyson equation, 
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has the following nontrivial solution: 
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Note that, when the value of coupling constant a lies in the region a < tt/6, X is real and the system is in 
the symmetric phase, while the value of a lies in the region a > tt/6, X is complex so that the solution of 
the Schwinger-Dyson equation is of the oscillatory function and the system is in the broken phase Q . The 
critical value ac in this case is tt/6, which is half of the well known value tt/3 in Landau gauge, because of 
the instantaneous exchange approximation and the simplification of the logarithm. 



The evaluation of the critica l val ue in the case including the full term of Eq.(A.2) is made as follows. 
The Schwinger-Dyson equation (A.l) with the expansion Eq.([A.2|) reads 
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0.475 



Figure 10: l/(ln(Af (p)/A))2 vs. a with a fitted line. The critical value is given by a^ = 0.44477 ± 0.00053. 



and we take the solution of the form M{p) ex p^, which is valid in the above simplified case. Note that the 
value of A is restricted by —2 < ReA < due to the condition that the integral in Eq.(A.(;) should be finite. 
Then the Schwinger-Dyson equation reduces to 



3a 



2(A+1) 
by performing the integral and using the formula 



^(A + 1) ^ 
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(A.8 



The real solution of Eq.(A.7) is easily found by considering the inclinations of tan7r(A + l)/2 and 2(A + l)/3 
at A = — 1. The result is that the nontrivial solution A is real when < a < 4/(37r) which corresponds to 
the symmetric phase, whereas Eq.( A.7 ) has complex solutions when a > 4/(37r) which corresponds to the 
broken phase. Hence the critical value of coupling constant in this case is 



4 
3^ 



(A.9) 



The exact value of critical coupling constant can be evaluated through numerical calculation, by fitting 
data with the function in the form of Miransky scaling law, and we find that the result is ac ~ 0.445. The 
numerical value 0.445 is in discrepancy with our theoretical prediction 4/(37r). It is easy to understand that 
this discrepancy is due to the fact that we made a linearization approximation. One should remember that 
the mass function in the denominator in Eq.(p3Q) is neglected in Eq.(A.l). This means that the exact value 
of ac must be slightly larger than 4/(37r) and would be in the range 



4 ^ TT 

37r '^'^ 6 
It is easy to see that our numerical result satisfies this speculation. 



(A.IO) 
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